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We microscopically study the effect of the magnetic field (Zeeman splitting) on the super- 
conducting state in a model for quasi-one-dimensional organic superconductors (TMTSF)2X. 
We investigate the competition between spin singlet and spin triplet pairings and the Fulde- 
Ferrell-Larkin-Ovchinnikov(FFLO) state by random phase approximation. While we studied 
the competition by comparison with the eigenvalue of the gap equation at a fixed temperature 
in our previous study (Phys. Rev. Lett. 102 (2009) 016403), here we obtain both the Tc for each 
pairing state and a phase diagram in the T(temperature)-/i2(field)-K; (strength of the charge 
fluctuation) space. The phase diagram shows that consecutive transitions from singlet pairing 
to the FFLO state and further to Sz ~ 1 triplet pairing can occur upon increasing the magnetic 
fleld when 2kF charge fluctuations coexist with 2A:_f spin fluctuations. In the FFLO state, the 
singlet d-wave and = triplet /-wave components are strongly mixed especially when the 
charge fluctuations are strong. 

KEYWORDS: organic conductor, quasi-one-dimensional system, (TI\/ITSF)2X, superconductivity, spin 
triplet pairing, FFLO state, Zeeman effect 



1. Introduction 

The superconducting state of quasi-one- 
dimensional (QID) organic conductors (TMTSF)2X 
(TMTSF=tetramethyl-tetraselenafulvalene, X=PF6, 
CIO4 etc.) has been an issue of great interest. From 
the discovery of the first organic superconductor 
(TMTSF)2PF6, various studies have been performed 
both experimentally and theoretically. ^"-'^'^ Previous 
studies for the NMR relaxation rate and 
the impurity effect^^"^^ have strongly suggested the 
possibility of anisotropic superconductivity where the 
nodes of the superconducting gap intersect the Fermi 
surface, although a thermal conductivity measurement 
has suggested the absence of nodes on the Fermi surface 
in (TMTSF)2C104.22 

Further experiments concerning the pairing symme- 
try have suggested the possibility that the pairing state 
in (TMTSF)2>'^ may be even more fascinating. The 
NMR Knight shift measurements for (TMTSF)2PF6 and 
(TMTSF)2C104 have shown that the Knight shift is un- 
changed across the superconducting critical temperature 
Tc. 23-25 The upper critical field Hc2 for (TMTSF)2PF6 
and (TMTSF)2C104 has been observed to exceed the 
Pauli paramagnetic limit iJp.^^-^s xhese experiments 
suggest the possibility of spin triplet pairing and/or the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, ^^-^o jj^ 
which the Cooper pairs formed as {k-\- |, — fc + J.) 
have a finite center of mass momentum Q^. 

Very recent experiments show more interesting re- 
sults. The NMR experiment for (TMTSF)2C104 shows 
that the Knight shift changes across Tc when the mag- 
netic field is small, but it is unchanged when a high 
magnetic field is applied. The Hc2 measurements for 
(TMTSF)2C104 show the possibility of two or three dif- 
ferent pairing states. In the intermediate field regime. 



superconductivity is easily destroyed by tilting the mag- 
netic field (between Hp and about 4T) from the con- 
ductive a-b plane, while in the high field regime, super- 
conductivity is sensitive to the broadening of the im- 
purity scattering potential of the nonmagnetic impurity, 
namely, if the broadening of the impurity potential is 
large, the upturn curve of the critical temperature van- 
ishes. 33 These experiments suggest that spin singlet 
pairing occurs in the field regime lower than the Pauli 
limit, but spin triplet pairing and/or the FFLO state 
occurs in the higher field regime. 

Theoretically, various studies on the pairing state in 
(TMTSF)2v'^ have elucidated not only the unconven- 
tional pairing state with a nodal gap function^^' 34-38 
but also the possibility of the spin triplet pairing3^-^'* 
and/or the FFLO state. ^'^"^^ In particular, we have pre- 
viously shown that the spin triplet "/-wave" pairing can 
compete with the spin singlet "d-wave" pairing in QID 
systems®" when 2kp spin fluctuations coexist with 2kp 
charge fluctuations since the Fermi surface is discon- 
nected in the fe-direction.3^-''i In fact, the coexistence 
of 2kp charge density wave and 2kp spin density wave 
in the insulating phase has been observed by the diffuse 
X-ray scattering experiments in (TMTSF)2PF6.*^^'®2 a 
similar conclusion concerning the pairing state competi- 
tion has been reached using the renormalization group 
technique. ^2. 43 niethod for identifying spin-triplet 

/-wave pairing, tunneling spectroscopy®3. 64 ^^le mid 
gap Andreev resonant state®^'®® and Josephson effect®^ 
have been proposed. In particular, the experiment of the 
proximity effect in the junctions with a diffusive normal 
metal is promising since an anomalous proximity effect 
with a zero-energy peak in the density of states, specific 
to spin-triplet superconductor junctions, has been pre- 

dicted.68-73 
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There have also been various studies for the pair- 
ing state in the magnetic field. The possibility of the 
FFLO state in finite magnetic field in a QID model for 
(TMTSF)2X has been suggested in several studies. ^^"^^ 
The possibility of field-induced spin triplet pairing has 
also been discussed by a phenomenological theory and a 
renormalization technique. ^^"^^ Recently, we have micro- 
scopically studied the magnetic field effect on the pairing 
state in QID systems. We found that the 5^ = 1 triplet 
pairing mediated by 2k p spin-|-2fci;' charge fluctuations is 
strongly enhanced by the magnetic field and showed the 
temperature-magnetic field phase diagram indicating the 
competition between the singlet and triplet pairings.'''' 
We further found that the spin singlet, triplet, and FFLO 
states are closely competing, and the Sz = Q triplet com- 
ponent is strongly mixed with the singlet component in 
the FFLO state. There, the pairing state competition 
has been studied by comparing the eigenvalue of the lin- 
earized gap equation in the space of Vy (strength of the 
charge fluctuation) and (magnetic field). 

The FFLO state has recently been studied actively not 
only in QID but also in general systems. ^^'^^ Previous 
theoretical studies have revealed various properties of the 
FFLO superconductivity from the viewpoint of (i) the or- 
bital effect, (ii) the impurity effect, and (iii) the 
anisotropy of the system.^'^"^'^^ One of the interesting as- 
pects of the FFLO state is parity mixing, i.e., even and 
odd parity pairings can be mixed to stabilize the FFLO 
state, which has been shown in phenomenological theo- 
ries.^°^'^°^ Recent microscopic studies have also shown 
that the = triplet pairing is mixed with singlet 
pairing in the FFLO state of the Hubbard model on the 
two-leg ladder-type lattice, the square lattice, -"^"^ 
and the QID extended Hubbard model. "''^ Yanase has 
pointed out that the parity mixing stabilizes the FFLO 
state, even in the vicinity of the quantum critical point, 
where the quasi-particle lifetime decreases owing to the 
scattering caused by spin fluctuations. In addition to 
these works, superconducting properties of the FFLO 
state have been studied theoretically. ^"^"^^^ 

Recent experiments strongly show the possibility 
of the FFLO state with an anisotropic gap function 
in CeCoIns.^^""^'^" Other candidate materials ex- 
hibiting the FFLO state are quasi-two-dimensional 
(Q2D) organic materials, such as A-(BETS)2-'^ 
(BETS=bisethylenedithio-tetraselenafulvalene, 
X=GaCl4,i3^ and FeCk^^"""*) and k-(BEDT- 
TTF)2Cu(NCS)2 (BEDT-TTF=bisethylenedithio- 
tetrathiafulvalene),^'*'''"^''^ and also a QID one 
(TMTSF)2C104.^^"^^ These materials have stimu- 
lated extensive studies in this field. The FFLO state 
attracts us not only in the field of superconductivity 
or superfiuidity in condensed matter but also in the 
quantum chromodynamics''^ and the ultracold fermionic 
atom gas.^'^**' ^^'^ 

Given the above background, in this study, we investi- 
gate the pairing competition between the spin singlet and 
spin triplet pairings, and the FFLO state of the supercon- 
ductivity mediated by spin and charge fluctuations in a 
QID extended Hubbard model for (TMTSF)2X by ran- 
dom phase approximation (RPA).^*°~^^^ While the com- 



petition was studied only by comparison with the eigen- 
value of the gap equation at a flxed temperature as in- 
dicated in ref. 54, here we calculate the superconducting 
temperature Tc for each pairing state. This enables us to 
obtain a phase diagram in the T(temperature)- ^.^(field)- 
Fj^ (strength of the charge fluctuation) space, where we 
find that (i) consecutive transitions from singlet pairing 
to the FFLO state and further to = 1 triplet pairing 
can occur upon increasing the magnetic field in the vicin- 
ity of the SDW-I-CDW phase, and (ii) the enhancement 
of the charge fluctuations leads to a significant increase 
in parity mixing in the FFLO state, where the 52 = 
triplet/singlet component ratio in the gap function can 
be close to unity. 

2. Formulation 

The extended Hubbard model for (TMTSF)2X [Fig. 
1(a)] that takes into account the Zeeman effect is given 
as 



+ 



E 

i,j.(T,(T' 



VijUi. 



(1) 



Here, Uj^ = tij +hzSgii{a)dij, where the hopping param- 
eters Uj considered are the intrachain (a-axis direction in 
(TMTSF)2X) nearest-neighbor tx and the interchain (6- 
axis direction) nearest- neighbor ty, = 1.0 being taken 
as the energy unit. U is the on-site interaction, and Vij 
are the off-site interactions: V^, Vx2, and V^s are the 
nearest-, next-nearest, and 3rd- nearest-neighbor interac- 
tions within the chains, and Vy is the interchain interac- 
tion. Note that we ignore the orbital effect, assuming that 
the magnetic field is applied parallel to the conductive 
x-y plane, assuming a sufficiently large Maki parameter. 
(Since we neglect the orbital effect, the direction of the 
magnetic field within the x-y plane is irrelevant within 
our approach.) 

The bare susceptibilities, consisting of bubble- type and 
ladder-type diagrams, are written as 



-i^ /(ea(fc + g))-/(Ug)) ^ (2) 



N 



^a{k + q) - ^a{q) ' 
^a{k + q) - ^a{q) ' 



(3) 



where ^cr{k) is the band dispersion that takes into ac- 
count the Zeeman effect measured from the chemical po- 
tential fi and /(^) is the Fermi distribution function. 

Within RPA that takes into account the magnetic field 
parallel to the spin quantization axis z,^^'^^ the longitu- 
dinal spin and charge susceptibilities are given by 

1 

2' 



Xs 



= Tiix^^ +x^^ -x^^ 



Xch = -^(x^^ + x^^ +x^^ + x 



(4) 
(5) 
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Fig. 1. (Color online) (a) Model adopted in this study, (b) 
Schematic figure of the gap for d-wave (left) and /-wave (right), 
where blue dashed lines indicate the nodes of the gap, and red 
solid curves indicate the disconnected Fermi surface. 



where 

X""{k) = [l + ^f{k)V{k)]xr{k)IA{k), (6) 

X^'ik) = -xr{k)[U + V{k)]xl''{k)IA{k), (7) 

A{k) = [l + x'',^{k)V{k)][l + xr{k)V{k)] 

-[U + V{k)fxZ''{k)xr{k). (8) 
The transverse spin susceptibiHty is given by 

Xt-{k) 



xTiik) = 



(9) 



l-Uxrik) 

where we ignore the ofF-site repulsions because it is dif- 
ficult to treat the effect of the off-site repulsions on the 
ladder diagrams. 

The pairing interactions from the bubble and ladder 
diagrams are given by 



^bu'L(^) = 



U + V(k) + ^xi;{k) 



[U + 2V{k)Y 



■Xch(fc), 



VilUk) = U\ti{k), 

V^:dk) = V{k)-2[U + V{k)]V{k)x'''{k) 
-V{kfx^%k) 
-[U + V{k)fx'''{kl 

v{L^{k) = 0. 



(10) 

(11) 



(12) 
(13) 



The linearized gap equation for Cooper pairs with the 
total momentum 2Qc {Qc represents the center of mass 



momentum) is given by 
1 

iV 



A5:'^--'(fc) = 1 Et^bub - 9) + Cd'(fc + q)] 



x ^^^f;»-/;-f-'(7-» ^-'(.), (14) 

where q± = q ± Qc, 'f'"' (fc) is the gap function, and 
X'q is the eigenvalue of this linearized gap equation. The 
center of mass momentum Q^, which gives the maximum 
value of A^*, lies in the x-direction,^'*'^^'!''^ while A^'^ 
takes its maximum at = (0, 0) because the electrons 
with the same spin can be paired as (fccr, —ka) for all k. 

We define the singlet and 5"^ = triplet components 
of the gap function in the opposite spin pairing channel 
as 



</'ss(fc) = 

(^STo(fc) = 



(pfi(fc) + (^J'T(fc) 



(15) 



In our calculation, the spin singlet and triplet compo- 
nents of the gap function in the FFLO state are essen- 
tially the d-wave and /-wave, respectively, as schemat- 
ically shown in Fig. 1(b); thus, we write the singlet 
(5*2 — triplet) component of the FFLO gap 9333 (v'ST") 
in eq. (15) as (/sssd ((PST/'Oi where SS(i(ST/°) stands for 
spin singlet d-wave (spin triplet /-wave with = 0) 
pairing. The eigenvalue of each pairing state is deter- 
mined as follows. Ag^ with ~ (0, 0) gives the eigen- 
value of the singlet d-wave pairing Xssd [Sz = triplet 
/-wave pairing Ast/o) 'PSTf = {(fissd = 0), while A^"^ 
with Qc ^ (0,0) gives Afflo- >^q[ with Q^ = (0,0) 
gives the eigenvalue for the spin triplet /-wave pairing 
with Sz = +1 {Sz = —1) AsT/+i (AsT/-i)- The above- 
mentioned results of the determination of the eigenvalues 
are listed in Table I. 



Table I. Results of the determination of the eigenvalue of the 
linearized gap equation Xq" . 





Center of mass momen- 
tum and paired spins 


Pairing symmetry or 
SC state 






singlet d-wave (Agg^j) 
for ipgTfO (fc) = 




= 0, cr ^ <t' 


Sz = triplet /-wave 
(AgT/o) for ^ssd (k) = 





= 0, cr = cr' 


5z = ±1 triplet /-wave 




^ 0, <7 7^ a' 


FFLO state (AppLo) 




^ 0, <7 = <7' 


not dominant state 



Although RPA is quantitatively insufficient for dis- 
cussing the absolute value of Tc, we expect this approach 
to be valid for studying the competition between differ- 
ent pairing symmetries. In this paper, we fix the hopping 



parameters as 



1.0 and 



electron interactions as, U — 1.7. 



0.2, and the electron- 
K = 0.9, 142 = 0.45, 



and Vx3, = 0.1, and vary Vy. Since the dimerization of 
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TMTSF molecules is very small in (TMTSF)2X com- 
pounds, we ignore the dimerization and fix the band 
filling as 71 — 1.5 (3/4 filling), where n = number of 
electrons/number of sites. 1024x128 fc-point meshes are 
taken, where we take a large number of meshes since 
the center of mass momentum Q^, which gives the max- 
imum value of the FFLO state, lies in the a;-direction. 

3. Results 

3. 1 Center of mass momentum and the gap function 

In this section, we study the nature of the FFLO state 
in our model. Let us first study the center of mass mo- 
mentum at which the FFLO state is most stabilized. 
The optimum that most stabilizes the FFLO state 
can be determined as at which the eigenvalue of the 
gap equation is maximized. In the following results, we 
set the interchain off-site interaction as Vy — 0.35. Aq^ 

with = (Qca;: Qcy) given in units of 7r/512 for 
the x-direction and 7r/64 for the y-direction. Figure 2 
shows the eigenvalue of the linearized gap equation in the 
opposite-spin pairing channel A^^ as a function of the x- 
component of the center of mass momentum Q^^ for var- 
ious Q^y. When the magnetic field is small {hz — 0.01), 
the pairing state with = (0, 0) dominates over other 
finite momentum states, as seen in Fig. 2(a). For a larger 
magnetic field {h^ = 0.03), a finite momentum pairing 
state with Q^^ = 3 and Q^y = dominates over other 
states in the opposite-spin pairing channel, as shown in 
Fig. 2(b). When the magnetic field is increased up to 
hz = 0.06, a finite momentum pairing state with Q^^ — 7 
and Q^y = is the most dominant, but the eigenvalue 
Ag^ itself decreases, as shown in Fig. 2(c). Further study- 
ing other hz cases, we find that the most dominant center 
of mass momentum lies in the x-direction,^^'^^' and 
the magnitude of the center of mass momentum increases 
with increasing magnetic field. 

The direction of the center of mass momentum vector 
can be understood from the Fermi surface split by the 
Zeeman effect shown in Fig. 3. The electron pair part, 
i.e., the particle-particle susceptibility, in the linearized 
gap equation is rewritten as 

f{Uq + Qc))-f{-(.'{-q + Qc)) 



i.a{q + Qc)+ia'{-q + Qc) 

iSn) Ga' {-q + Qc 



(16) 



where (,(j{—k+Qc) is the same as £,a{k—Qc) since ^(j(fc) ~ 
^ai—k) is satisfied. If the a spin electron energy at the 
wave vector q + Qc and the a' spin electron energy at 
the wave vector —q + Qc are close to the Fermi energy, 
the denominator is small and eq. (16) can take a large 
value. For a quasi-one-dimensional system, the number 
of wave vectors q that satisfies such a condition becomes 
the largest when the vector Qc is in the k^ direction 

Next, we study the gap functions normalized by the 
maximum value of the singlet component gap function 
in the FFLO state. We set the parameters as hz — 0.03, 
Vy = 0.35, and T = 0.012, where the FFLO state which 
has the finite center of mass momentum as {Qcx, Qcy) = 
(3, 0) is the most dominant, as described later. Note that 
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Fig. 2. (Color online) Qj.^-dependence of the eigenvalue in the 
opposite-spin pairing channel, A^"^, for (a) hz = 0.01, (b) hz = 
0.03, and (c) hz = 0.06 at Vy = 0.35. 



the 5^ = ±1 triplet pairings always have the maximum 
value of the eigenvalue Xq' at = (0, 0), as mentioned 
previously. As shown in Figs. 4(a) and 4(b), the singlet 
component of the gap function in the FFLO state is the 
d-wave and the Sz = triplet component is the /-wave. 
The maximum value of the Sz — triplet gap compo- 
nent in the FFLO state almost reaches unity. Thus, the 
singlet c?-wave component and the Sz = triplet /-wave 
component strongly mix in this FFLO state. The gap 
function in the = ±1 triplet pairings has the /-wave 
form shown in Figs. 4(c) and 4(d). 
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kx 

Fig. 3. (Color online) The small purple arrow in the fcx-direction 
denoting QppLO schematically shows the direction of the center 
of mass momentum vector Q^. in the FFLO state. The black thin 
solid cuarves schemetically represent the Fermi surface in zero 
field. The red thick solid (green thick dashed) curves schemati- 
cally represent the Fermi surface split by the Zeeman effect. The 
red and green filled circles are the particle on each Fermi surface 
in the presence of the field, and the gray filled squares are the 
particle in zero field. 




Fig. 4. (Color online) Gap function for (a) singlet component, (b) 
Sz = triplet component in the FFLO state with Qcx = 3 and 
Qcy = 0, (c) Sz = +1 triplet state, and (d) Sz = —1 triplet 
state, where black solid curves represent the Fermi surface, and 
the green dashed curves represent the nodes of the gap. The 
parameters are hz = 0.03, Vy = 0.35, and T = 0.012. 

The appearance of the c?-wave gap in the singlet com- 
ponent and the /-wave gap in the 5*2 = triplet com- 
ponent in the FFLO state is understood as follows. In 
zero field, the singlet c?-wave pairing mediated by the 
2kF spin fluctuations is favored in the QID Hubbard 
model, namely, the large pairing interaction due to the 
2kp spin fluctuations stabilizes the spin singlet c?-wave 
pairing. '^''"^^'^^ Moreover, the coexistence of 2kp charge 
fluctuations, which is induced by the second- nearest- 



neighboring repulsive interaction, favors the triplet /- 
wave pairing in the QID extended Hubbard model at 
quarter filling. '^^^'^^ The reason why the spin triplet /- 
wave pairing can compete with the spin singlet c?-wave 
pairing in the QID extended Hubbard model is (i) the 
contribution of the 2kp charge fluctuations in the pairing 
interaction enhances the spin triplet /-wave pairing and 
suppresses the spin singlet d-wave pairing, and (ii) / and 
d-wave pairings have the same number of gap nodes in- 
tersecting the Fermi surface due to the disconnectivity of 
the Fermi surface (quasi-one-dimensionality). The above 
mechanism is valid even in the presence of the magnetic 
field, but more importantly, the spin triplet /-wave pair- 
ing mediated by the 2kp spin -f 2kp charge fluctuations 
can be enhanced by applying the magnetic field since the 
bubble-type diagram enhanced by the field contributes 
to the pairing interaction without being paired with the 
bubble-type diagram, which is suppressed by the field. ^'^ 
Actually, our previous work shows a clear correlation be- 
tween the Sz = triplet ratio in the FFLO state and 
the ratio of the eigenvalue between the 5'z = triplet 
and singlet pairings obtained by the formulation of sep- 
arating the singlet and Sz = Q triplet channels.^'* From 
the above, we can understand not only the appearance 
of the c?-wave (/-wave) gap in the singlet {Sz = triplet) 
component of the opposite-spin pairing channel and the 
/-wave gap in the parallel-spin pairing channel, but also 
the large parity mixing of the singlet and Sz ^ ^ triplet 
components in the FFLO state. 

Figure 5 shows the parity mixing ^sTf°lvsSd in the 
opposite-spin pairing channel as a function of the x- 
component of the center of mass momentum Qcx for sev- 
eral Qcy Note that we need to bear the dependence 
of the eigenvalue A^^, as shown in Fig. 2, in order to 
see the dependence of the parity mixing rate be- 
cause the most dominant pairing state in the opposite- 
spin pairing state is determined by the value of Aq^ . For 
instance, we have seen in Fig. 2(a) that the singlet d- 
wave pairing, i.e., the opposite-spin pairing state with 
iQcxTQcy) = (0,0), is the most dominant in the small 
magnetic field regime. In Fig. 5(a), the parity mixing rate 
for hz = 0.01 ipsTf/'Pssd is zero at (Qcx, Qcy) = (0,0). 
Therefore, no 5*2 = triplet /-wave component is present 
in this pairing state, and the opposite-spin pairing chan- 
nel is a purely spin singlet d-wave. For hz — 0.03, we have 
seen in Fig. 2(b) that the FFLO state with Qcx — 3 and 
Qcy = is dominant. As shown in Fig. 5(b), the parity 
mixing rate for Qcx = 3 and Qcy — takes a large value 
ifSTf/^ssd — 0.8. For hz — 0.06, where the FFLO state 
with Qcx = 7 and Qcy = is dominant (Fig. 2(c)), the 
parity mixing rate increases, i.e., fsTf/vsSd — 1-0, as 
shown in Fig. 5(c), which means that the singlet d-wave 
component and the Sz ~ triplet /-wave component are 
strongly mixed in this FFLO state (provided this state is 
actually realized). The strong parity mixing in the FFLO 
state can be understood as a consequence of the breaking 
of the spacial inversion symmetry in the superconduct- 
ing state. Previous theoretical studies have shown that 
the parity mixing with the singlet and triplet pairings 
stabilizes the FFLO state more when only the singlet 
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Fig. 5. (Color online) -dependence of the parity mixing in the 
opposite-spin pairing channel, "PsT/o/VSSdi for (^) f'-z = 0.01, 
(b) = 0.03, and (c) = 0.06 at Vy = 0.35. 



component is considered. 

3.2 Temperature dependence 

Next, we investigate the temperature dependence of 
the eigenvalue, Xq' , in both the opposite- and parallel- 
spin pairing states. We have confirmed that the cen- 
ter of mass momentum that most stabilizes the FFLO 
state is unchanged upon lowering the temperature for 
a fixed magnetic field. For = 0.01, the eigenvalue 
Assd = -^0^=0 ^P^^ singlet d-wave pairing reaches 

unity, as shown in Fig. 6(a). In this small magnetic field 



Fig. 6. 

\<7a 



(Color online) Eigenvalue of the linearized gap equation, 
plotted as a function of the temperature T for (a) hz = 



SSd and ST/±l have = 0, and the FFLO state has a finite 
that maximizes the eigenvalue of the opposite-spin channel 
in Fig. 2. 



regime, the FFLO state is absent, as shown in Fig. 2(a). 
For hz — 0.03, the singlet c?-wave pairing is suppressed 
and the eigenvalue Afflo = -^q^/o '^^ ^^'^ FFLO state 
with Qcx — 3 and Qcy = reaches unity as seen in 
Fig. 6(b). For h,, = 0.06, the FFLO state with Qcx = 7 
and Qcy = does not develop much upon lowering the 
temperature, while the eigenvalue Asx/+i = ^cf =o 
the 5*2 = 1 triplet /-wave state reaches unity, as shown 
in Fig. 6(c). The eigenvalue of the singlet c?-wave and 
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5*2 = —1 triplet /-wave pairings remains small even in 
the low temperature regime. 

3. 3 Calculated phase diagram 

We now obtain a phase diagram in the temperature 
T versus the magnetic field space for several values 
of interchain off-site interaction Vy (which controls the 
strength of the charge fluctuations). Figure 7(a) shows a 
plot of the critical temperature Tc against the magnetic 
field hz for Vy — 0.35, where the 2kp charge fluctuations 
are slightly weaker than the 2kp spin fluctuations. The 
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Fig. 7. (Color online) (a) Calculated phase diagram in hz-T space 
for Vy = 0.35, where the green dashed line indicates the Tc for 
the spin singlet d-wave, the red solid line represents that for 
the FFLO state, and the blue dotted line indicates that for the 
Sz = 1 spin triplet /-wave. The spin singlet d-wave is omitted 
as SSrf and the 5z = 1 spin triplet /-wave is ST/^'^. The same 
notation is used in Figs. 8 and 9. (b) Schematic figure of the or- 
bital pair breaking effect on the superconducting phase diagram 
in T-hz space, where black solid arrows schematically represent 
the orbital pair breaking effect. 



critical temperature in zero field is Tc ~ 0.012 and the 
estimated value of Pauli's paramagnetic field is ~ 
0.03. We see that a consecutive transition from singlet 
pairing to the FFLO state and further to 5*2 = 1 triplet 
pairing occurs upon increasing the magnetic field. 

This consecutive pairing transition can be understood 
as follows. It is known that the FFLO state can be stabi- 
lized by the quasi-onc-dimensionality, namely, the nest- 



ing of the Fermi surface.^'*' Thus, the quasi-one- 
dimensionality of the present model is one of the origins 
of the transition from the d-wave to the FFLO state. The 
origin of the pairing transition from the FFLO state to 
the Sz = ^ triplet pairing is understood by our previous 
study, where we have shown that the triplet pairing due 
to the coexisting 2k p spin and 2fci? charge fluctuations 
is strongly enhanced by the direct contribution of the 
unpaired bubble diagram enhanced by the field. 

Here, we emphasize that we ignore the orbital pair 
breaking effect in this study because our aim in this work 
is to study the competition between the singlet, FFLO, 
and triplet pairings in the case when the magnetic field 
is applied in the conductive plane, i.e., the a-b plane of 
(TMTSF)2v'C. For discussing the above pairing competi- 
tion, the Zeeman splitting effect is essential for the FFLO 
state; thus, we ignore the orbital pair breaking effect at 
the beginning. Although this effect is small in applying 
the magnetic field parallel to the conductive plane, the 
orbital pair breaking effect is present in actual materi- 
als. Furthermore, previous studies have shown that the 
orbital pair breaking effect is important in discussing the 
FFLO superconductivity.''^"*'' 

As shown in Fig. 7(a), there seems to be a reentrance 
from the superconducting state to another superconduct- 
ing state intervened by the normal state. However, it is 
much more reasonable to consider that this reentrance 
does not actually occur owing the presence of the or- 
bital pair breaking effect. If this effect is taken into ac- 
count, not only the FFLO state but also the singlet 
and triplet pairing states should strongly be suppressed 
upon increasing the magnetic field. Figure 7(b) shows a 
schematic figure of the effect of the orbital pair break- 
ing, where the Tc obtained (without the orbital effect) 
in Fig. 7(a) (thin curve) is suppressed down to the thick 
curve. The thick curve in Fig. 7(b) is reminiscent of the 
experimental T-H phase diagram^^"^*''^^''^^ in that the 
Tc curve makes an upturn from nearly above the Pauli 
limit. 

Next, we study the effect of the interchain interac- 
tion Vy on the phase diagram in the temperature T 
versus the magnetic field hz space. Figure 8(a) shows 
the critical temperature Tc at each magnetic field hz for 
Vy = 0.38. The magnetic field at which the transition 
from the FFLO state to the Sz — ^ triplet /-wave pairing 
occurs is smaller than that in the Vy = 0.35 case. The 
critical temperature Tc for Vy = 0.32 is shown in Fig. 
8(b), which shows that the magnetic field at which the 
FFLO state gives way to the Sz — I triplet /-wave pair- 
ing is larger than those in the previous phase diagrams 
shown in Figs. 7(a) and 8(a). The difference between the 
two phase diagrams is due to the fact that the 2kF charge 
fluctuations enhance the triplet /-wave pairing; thus, the 
FFLO state appears only in a small parameter regime in 
between the d- and /-wave pairings. 

Summarizing the above-mentioned features, we show 
the phase diagram in T-Vy-hz space in Fig. 9. When Vy 
is small and thus the 2k p spin fluctuations are domi- 
nant over the 2kp charge fluctuations, Tc decreases and 
the transition from the spin singlet d-wave to the FFLO 
state occurs upon increasing hz- In this FFLO state, the 
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Fig. 8. (Color online) Calculated phase diagram in h^-T space 
for (a) Vy = 0.38 and (b) Vy = 0.32, where the notation is the 
same as that in Fig. 7(a). 




Fig. 9. (Color online) The critical temperature is shown in the 
Vy-hz plane. The value of Tc is plotted in the vertical axis and 
represented by contours. Green dashed lines represent the critical 
temperature for the singlet d-wave pairing, red solid lines are the 
Tc for the FFLO state, and the blue dotted lines are the Tc for 
the Sz = 1 spin triplet /-wave pairing, respectively. 



strong parity mixing with the spin singlet c?-wave com- 
ponent and the Sz = spin triplet /-wave component 
occurs. In the large Vy regime, the 2kF charge fluctua- 
tions compete with the 2kF spin fluctuations, and the 
consecutive pairing state transition from the spin singlet 
d-wave to the FFLO state and further to the Sz = I spin 
triplet /-wave upon increasing the hz at the critical tem- 



perature Tc occurs. The Tc enhancement of the Sz — I 
spin triplet /-wave pairing in the large hz regime can be 
understood by our previous work.^"^ 

4. Conclusion 

We have studied the competition between spin sin- 
glet, triplet, and FFLO superconductivities in a model 
for (TMTSF)2X by applying the RPA method and solv- 
ing the linearized gap equation within the weak coupling 
theory. We find the following: 

(i) consecutive pairing transitions from singlet pairing 
to the FFLO state and further to Sz = 1 triplet pair- 
ing can occur upon increasing the magnetic field in the 
vicinity of the SDW+CDW coexisting phase. 

(ii) in the FFLO state, the Sz ~ spin triplet pairing 
component is mixed with the spin singlet pairing com- 
ponent, thus resulting in a large parity mixing. 

Recent experiments for (TMTSF)2C104 suggest differ- 
ences in superconducting properties in the low and high 
field regimes. The Knight shift study shows the pres- 
ence of low field and high field pairing states, where the 
former is the spin singlet pairing and the latter is the 
FFLO state or the spin triplet pairing. '^^ The upper crit- 
ical field studies have shown that only the clean sample, 
or more strictly, samples where the broadening of the 
nonmagnetic impurity is small, exhibits an upturn of the 
critical temperature curve in the high field parallel to 
the a axis regime above 4T; thus, the high field super- 
conducting state is sensitive to the impurity content or 
the anisotropy of the impurity scattering potential. '^^ Be- 
tween 4T and the Pauli limit around 2.5T, there seems 
to be a different high field pairing state, in which su- 
perconductivity is stable against the impurities, but it 
is very sensitive to the tilt of the magnetic field out of 
the a-b plane. The bottom line of these experiments is 
that there may be three kinds of pairing states, i.e., one 
low field state and two high field states. The correspon- 
dence between these experimental observations and the 
present study is not clear at the present stage, but the 
appearance of the three kinds of pairing states is indeed 
intriguing. It would be interesting to further investigate 
experimentally the possibility and nature of two kinds of 
high field pairing states. 

One point that should be mentioned for 
(TMTSF)2C104 in particular is the presence of the 
anion ordering with the modulation wave vector 
Qao — which takes place near Tao — 24K 

when slowly cooled. Recent studies show that the anion 
ordering potential (Vao) is around 0.02t^}'^'^-^'^^ The 
anion ordering leads to a folding of the Brillouin zone 
in the fcy-direction (6-direction), and in that case, 
the d-wave (and also /-wave in the same sense) gap 
can become nodeless because the folded Fermi surface 
becomes disconnected near the nodes of the gap, as 
has been suggested by Shimahara.^^^ This effect is 
neglected in our present study, and its effect on the 
pairing symmetry competition is an interesting future 
problem. 

Another point to be mentioned is that in the present 
study, we do not take account of the retardation effect. 
By taking account of this effect, i.e., the frequency de- 
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pendence of the gap function, we can discuss the odd- 
frequency pairing state. It has been shown that 
odd-frequency pairing can be reahzed in a certain quasi- 
one-dimensional lattice. In particular, in the presence 
of non-uniformity, the odd-frequency pairing amplitude 
is ubiquitously generated. It is a future interesting 
problem to study the possible existence of odd-frequency 
pairing in quasi-one-dimensional organic superconduc- 
tors. 
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